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Abstract. ANTARES is a submarine neutrino tele- 
scope deployed in the Mediterranean Sea, at a depth 
of about 2500 m. It consists of a three-dimensional 
array of photomultiplier tubes that can detect the 
Cherenkov light induced by charged particles pro- 
duced in the interactions of neutrinos with the 
surrounding medium. Down-going muons produced 
in atmospheric showers are a physical background to 
the neutrino detection, and are being studied. In this 
paper the measurement of the Depth Intensity Re- 
lation (DIR) of atmospheric muon flux is presented. 
The data collected in June and July 2007, when the 
ANTARES detector was in its 5-line configuration, 
are used in the analysis. The corresponding livetime 
is 724 h. A deconvolution method based on a Bayesian 
approach was developed, which takes into account 
detector and reconstruction inefficiencies. Compari- 
son with other experimental results and Monte Carlo 
expectations are presented and discussed. 

Keywords: Cherenkov neutrino telescope, underwa- 
ter muon flux. 

I. Introduction 

The largest event source in neutrino telescopes is 
atmospheric muons, particles created mainly by the 
decay of it and K mesons originating in the interaction 
of cosmic rays with atmospheric nuclei. Although 
ANTARES ffl, El, 0, El, "looks downwards" in 
order to be less sensitive to signals due to downward 
going atmospheric muons, these represent the most 
abundant signal due to their high flux. They can be a 
background source because they can be occasionally 
wrongly reconstructed as upward going particles 
mimicking muons from neutrino interactions. On the 
other hand they can be used to understand the detector 
response and possible systematic effects. In this scenario 
the knowledge of the underwater // intensity is very 
important for any Cherenkov neutrino telescope and the 
future projects 0, Q. Moreover, it would also provide 
information on the primary cosmic ray flux and on the 
interaction models. 

II. Data and simulation samples 

The considered data sample is a selection of June 
and July 2007 data: only runs with good background 
condition^] are considered in the analysis. The livetime 

'Good run: averaged baseline rate below 120 kHz, burst fraction 
(due to biological activity) below 20%, muon trigger rate more than 
0.01 Hz and less than 10 Hz. 



MUPAGE generation parameters 




Min 


Max 


Shower Multiplicity 


1 


100 


Shower Energy (TeV) 


0.02 


500 


Zenith angle (degrees) 


95 


180 



TABLE I: Generation parameters set in the MUPAGE 
simulation. 



of the real data sample corresponds to 724 h. 

Atmospheric muons were simulated for the 5-line 
ANTARES detector. The equivalent livetime corre- 
sponds to 687.5 h. The Monte Carlo programs used in 
the analysis are the following: 

. Physics generator: MUPAGE program ©, g). It 
generates the muon kinematics on the surface of an 
hypothetical cylinder (can) surrounding the detector 
instrumented volume (see Tab. H). 

• Tracking and Cherenkov light generation: KM3 
program IfTOl . 

A dedicated program inserts the background in the 
simulation taking it from a real run. The Monte Carlo 
data are then processed by the trigger software, which 
requires the same trigger conditions as in the real data. 

Physical information is inferred from the triggered 
events (both Monte Carlo and real data) by a chi-square 
based track reconstruction program HI 111 . Each event is 
reconstructed as a single muon, even if it is a muon 
bundle. 

III. Cut selections based on the track 

RECONSTRUTION ALGORITHM 

A chi-square reconstruction strategy IfTTII is used in the 
analysis. Different fits, based on a chi-square minimiza- 
tion approach, are applied by the tracking algorithm: 

- a linear rough fit whose extracted parameters are used 
as starting point for the next refined fits; 

- a track fit which looks for a muon track; 

- a bright point fit which looks for a point light source 
as for example electromagnetic showers originated by 
muon interactions with matter. 

Particular interest in the analysis is given to the follow- 
ing quality parameters: 

• nline: number of lines containing hits used in the 
track fit algorithm; 

• nhit: number of hits (single or merged) used in the 
track fit algorithm; 

• Xt : normalized chisquare of the track fit. The 
smaller is its value the larger is the probability that 
the reconstructed track belongs to a muon event; 
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Efficiency(%) 
Real data 


Efficiency(%) 
MC data 


Purity(%) 
MC data 


No cut 


100 


100 


62 


Reconstructed track 


99 


99 


63 


nhit > 5* 


89 


94 


64 


Xt <3 


51 


54 


77 


xi>2 


50 


53 


78 



TABLE II: Efficiencies and Purities. The cuts are per- 
formed in sequence. 
*nhit > 5 applied only on SL events. 



• xt '■ normalized chisquare of the bright point fit. The 
larger is its value the larger is the probability that 
the reconstructed track belongs to an electromag- 
netic shower and not to a muon. 

Some cuts, based on these quality parameters, are 
necessary to improve the purity of the data sample. 

The hits used in the track reconstruction may belong 
only to one line (Single Line-SL event) or to more than 
one line (Multiple Line-ML event). The events detected 
with a single line usually have a well reconstructed 
zenith angle but undefined azimuth angle (if the line is 
perfectly vertical, the hit informations are independent 
from the azimuth angle of the track). The measurement 
of the Depth Intensity Relation is not strictly related with 
the azimuth angle and for this reason single line events 
are also considered here. 

The cuts are performed in sequence on the quality pa- 
rameters of the reconstruction program. The Efficiencj@ 
and the Purit}0 of the selected data set after each cut are 
presented in Tab. [TT] The first cut is needed in order to 
remove the events for which the reconstruction algorithm 
does not converge toward a definite value of the fitting 
parameters. 

IV. Depth Intensity Relation 

In the present section the quantities are given as func- 
tions of the zenith angle $ obtained from the unfolded 
real events. 

One method to derive the DIR is to compute the muon 
flux Ih (cos •&) as a function of the zenith angle $ at a 
fixed vertical depth ho in the sea. Once this distribution 
is known, it can be transformed into the DIR using the 
relation fl2l. Ifl3l: 



Iv(h) 



2/j (cos I?) • I cosi?| • k(cosi?) (1) 



cm 2 ■ sr 1 



where the subscript "V" stands for "Vertical events'" (i.e. 
cost? = —1) and h = /in./ cost? represents the slant 
depth (i.e. the distance covered in the sea water by 
muons, to reach the vertical depth ho with zenith angle 
direction ff). In the following ho = 2000 m corresponds 

2 The Efficiency is defined as the fraction of events surviving the 
cuts over all the reconstructed events. 

3 The Purity is defined as the fraction of events with a zenith 
reconstruction error A8 = \6t — m \ < 5° over the selected events. 
Applicable only to MC. 




-0.3 -0.2 -0.1 

COS0 t 

Fig. 1: Rh (cos8 t ). From Monte Carlo. Only statisti- 
cal errors are shown. 



to the sea depth of the highest ANTARES Photomul- 
tiplier Tubes (PMTs). The equation Q] is referred to 
as "flux verticalization": it transforms the muon flux 
7fc, (cos$) as a function of the zenith angle t? at the 
fixed sea depth ho, into the DIR Iv(h), i.e. the flux 
of the vertical muons as a function of the sea depth h. 
The | cost?| and the k(cosi?) factors are needed in order 
to take into account the zenith angle dependence of the 
atmospheric muon flux at sea level 

m, ma. 

The measured zenith distribution N R (cos9 m ) is ob- 
tained from the track reconstruction of the selected real 
events. The deconvolution procedure is a method to 
derive a true distribution from a measured one. In this 
work the goal is to transform the measured real data dis- 
tribution N^(cos 6 m ) into its parent angular distribution 
N R (cos t?), which represents the real events crossing the 
can surface during the considered experimental time: 

N R (cos9 rn ) — > (deconvolution) — ► N R (cos , &) 

(2) 

This is possible using the Monte Carlo simulations of 
the detector response. 

Several methods to unfold data exist. The approach 
that has been chosen consists in an iterative method 
based on Bayes' theorem proposed in lfl4l . 

Once the distribution N R (cos$) is known, it is possi- 
ble to derive the atmospheric muon flux Ih (cos $) at the 
fixed depth ho. From relation Q] the DIR can be finally 
written as in the following equation: 

N R (cosd) -rnhgicos-d) ■ R ho (cos ■&) 



Iv(h) 



AT ■ Af2 • A c (cost?) 
cost?) ■ k(cos^) [s~ 1 ■ cn 



(3) 



sr 



where the quantities in the equation are the followings: 
- AT = 2.61 • 10 6 s is the livetime of the considered 
real data sample. 

4 The sea level flux has a zenith angle dependence oc 1/ (cos i9 ■ 
re(cos$)) where the corrective factor is needed to take into consider- 
ation the Earth curvature. re(cosi?) can be considered equal to 1 for 
cosi? < -0.5. 
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Fig. 2: Average muon event multiplicity rnh (cos9t) at 
the fixed sea depth h = 2000 m for > 20 GeV. 
From Monte Carlo. Only statistical errors are shown. 
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Fig. 3: Relative difference between the DIRs obtained 
with the defined quality cuts and without quality cuts 
(see eq. 0. 



- AO = 2tt ■ 0.05 sr is the solid angle subtended by two 
adjacent zenith angle bins as considered in the analysis. 

- A c (cos$) is the generation can area as seen under the 
zenith angle $ (projection of a cylinder): 

A c (costf) = ttRI • I cost?| +2R c -H c -\ sintf| (4) 

R c = 511m and H c = 585 to are the radius and the 
height of the generation can. 

- N R {cosd) represents the number of muon events 
reaching the generation can surface during the consid- 
ered experimental time AT. 

- i?/, (cos $) is a correction factor needed to get the 
event flux at the sea depth ho from the event flux 
averaged on the whole can area. Rh (cos computed 
from Monte Carlo ($ = 8 t , where 9 t is the -"true"- 
generated zenith angle of the Monte Carlo muon event), 
is shown in Figure Q] 

- mh (cos $) is the average muon bundle multiplicity 
at the fixed sea depth ho = 2000 to. This quantity, 
computed from Monte Carlo ($ = 9 t ), is shown in 
Figure [2] 

- k(cos$) • cos($) are the correction factors fl2l . |[T3l 
introduced in eq. Q] 

V. Estimation of systematic uncertainties 

During MC simulation several input parameters are 
required to define the environmental and geometrical 
characteristics of the detector. Some of them are con- 
sidered as sources of systematic uncertainties. In ITT31 
the effect of the variation of the following quantities on 
the muon reconstructed track rate is considered. 

• Modifying by ±10% the reference values of the sea 
water absorption length, an almost negligible effect 
on the shape of the zenith distributions was noticed, 
while the absolute flux changed by +18%/ - 20%. 

• Decreasing and increasing the effective area of 
the ANTARES optical module (OM) by 10% with 
respect to the values used in the analysis, a change 
of about ±20% was observed in the muon flux. 



> The effect of the maximum angle between the OM 
axis and the Cherenkov photon direction allowing 
light collection was considered. Moving the cut- 
off of this OM angular acceptance the rate of 
reconstructed tracks change of about +35%/— 30%. 
Summing in quadrature the different contributions, a 
global systematic effect of about +45%/ — 40% can 
be considered as an estimate of the errors produced by 
uncertainties on environmental and geometrical param- 
eters. 

The obtained results also depend on the quality cuts 
performed on the data set. The unfolding algorithm 
depends on the relative ratio of Monte Carlo and real 
reconstructed events used in the analysis. As seen in 
Tab. the selections applied to the events have slightly 
different effects on the two data sets. In order to take 
into account this effect the unfolded DIR Iy(h) has been 
obtained without considering any cut but the first one 
which eliminates not reconstructed tracks. The relative 
difference K{h) between the two final DIRs is defined 
in the following equation: 

i;.(h)-I v (h) 



K{h) 



Iv(h) 



(5) 



The quantity, considered as a systematic uncertainty, 
depends on the slant depth and is shown in Figure[3] This 
uncertainty is summed in quadrature with the systematic 
error estimated above. 

In Figure g] the muon flux I ho (cos 1?) (E^ > 20 GeV) 
at 2000 m depth is shown. The Monte Carlo simulation 
from MUPAGE is also displayed. The result, as the 
Monte Carlo simulation, takes into account only muons 
with energy higher than 20 GeV because muons with 
lower energy are not able to trigger the detector. 

In Figure [5] the DIR Iv{h) is shown together with 
other experimental results. The Sinegovskaya parameter- 
ization (E^ > 20 GeV) [IT) and the Monte Carlo sim- 
ulation from MUPAGE are also shown. The ANTARES 
results are in reasonable agreement with the ones of the 
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Fig. 4: PRELIMINARY. Flux of atmospheric muons 
for Ef, > 20 GeV at 2000 m of sea depth (I h „ (cos J?)) 
with systematic uncertainties (the statistical uncertainties 
are negligible). The MUPAGE simulation is superim- 
posed. 



deep telescope prototypes DUMAND 03 and NESTOR 
Ifl9l and with the data of the Baikal EH), ED and 
AMANDA E2) collaborations. 

VI. Conclusions 

The aim of the presented analysis is the measurement 
of the muon flux at the depth of ANTARES and the 
derivation of the vertical component of the atmospheric 
muon flux as a function of the sea depth. The goal is 
also to assess the performance of ANTARES in detecting 
muons. The analysis has been performed on a selection 
of the experimental data of June and July 2007 when 
the ANTARES detector was in its 5-line configuration. 

Several quality cuts have been applied on the re- 
constructed events in order to improve their purity, in 
particular concerning the zenith angle reconstruction. 

An unfolding algorithm, based on an iterative method, 
has been applied on the selected experimental data in 
order to retrieve back the flux of atmospheric muons 
with Efj, > 20 GeV at the fixed sea depth ho = 2000 m. 
The experimental DIR was finally obtained. 

The results are in good agreement, within the uncer- 
tainties, with the experimental fluxes obtained by other 
Cherenkov telescopes. 
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